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ABSTRACT 

In this paper, we investigate the applicability of the astrometric method to 
the determination of the lens parameters for gravitational microlensing events 
toward both the LMC and the Galactic bulge. For this analysis, we investigate 
the dependency of the astrometrically determined angular Einstein ring radius, 
^(^e/^e,o); on fhe lens parameters by testing various types of events. In 
addition, by computing A(6*e/^e,o) for events with lensing parameters which 
are the most probable for a given lens mass under the standard models of 
Galactic matter density and velocity distributions, we determine the expected 
distribution of the uncertainties as a function of lens mass. From this study, 
we hnd that the values of the angular Einstein ring radius are expected 
to be measured with uncertainties A(6*e/6*e,o) ^ 10% up to a lens mass of 
M ~ 0.1 Mq for both Galactic disk-bulge and halo-LMG events with a moderate 
observational strategy. The uncertainties are relatively large for Galactic 
bulge-bulge self-lensing events, A(6'e/6'e,o) ~ 25% for M ~ 0.1 Mq, but they can 
be substantially reduced by adopting more aggressive observational strategies. 
We also hnd that although astrometric observations can be performed for most 
photometrically detected Galactic bulge events, a signihcant fraction (~ 45%) 
of LMG events cannot be astrometrically observed due to the faintness of their 
source stars. 
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1. Introduction 


Current microlensing experiments (Alcock et al. 1997a, 1997b; Ansari et al. 1996; 
Udalski et al. 1997; Alard & Guibert 1997) are detecting massive astrophysical compact 
objects (MACHOs) by monitoring the light variations of stars undergoing gravitational 
microlensing. The light curve of a microlensing event is related to the lensing parameters by 
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where A is the magnification, u is the lens-source separation, (3 is the impact parameter, 
tp is the time of maximum amplification, and is the Einstein ring radius crossing time 
(Einstein timescale). Of these lensing parameters, only the timescale provides information 
about the lens because it is the only one directly related to the physical parameters of the 
lens (lens parameters) by 
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where v is the lens-source transverse speed, te is the physical size of the Einstein ring 
radius, M is the mass of the lens, and Doi, Dig, and Dos are the separations between the 
observer, lens, and source star. However, due to the degeneracy in the lens parameters, the 
nature of the lens population is still very uncertain. 

The uncertainties in the lens parameters can be signihcantly reduced if one can measure 
the lens proper motion, p = 6 *e//e, where 9^ = r-^/Doi is the angular Einstein ring radius. 
The lens proper motions can be measured both photometrically (Gould 1994; Nemiroff & 
Wickramasinghe 1994; Witt & Mao 1994) and spectroscopically (Maoz & Gould 1994; Loeb 
& Sasselov 1995). However, measuring the proper motion with either of these techniques is 
only possible for several special classes of lensing events. 

The lens proper motions can be measured in a general way if the separation between 
the two source star images can be measured. As this image separation is generally on order 
1 pas, it can not be measured with current instrumentation. However, the proposed Space 
Interferometry Mission (Allen, Shao, & Peterson 1997, hereafter SIM) will have high enough 
positional accuracy to measure the astrometric displacement in the light centroid caused by 
gravitational microlensing (Walker 1995; Hpg, Novikov & Polnarev 1995; Paczyhski 1998; 
Boden, Shao, & Van Buren 1998). The astrometric shift of the source star image centroid 
is related to the lensing parameters by 


where T = [t — tp)/tE and the x and y axes represent the directions which are parallel and 
normal to the lens-source transverse motion, respectively. The trajectory of the apparent 



4 


source star image traces out an ellipse during the event (astrometric ellipse, Jeong, Han, 

& Park 1998). With the measured astrometric ellipse, one can determine the impact 
parameter of the event because the shape (i.e. the axis ratio) of the ellipse is related to 
the impact parameter 13. The importance of measuring the astrometric centroid shifts is 
that one can then determine the angular Einstein ring radius as the size of the astrometric 
ellipse, i.e., the semimajor axis, is directly proportional to 6-e. Once the angular Einstein 
ring radius is determined from the astrometric shifts, one can determine the lens proper 
motion from 6^ combined with the event timescale, which is determined from photometric 
monitoring of the event. 

Since astrometric centroid shifts of microlensing events exhibit various shapes depending 
on the combination of the lensing parameters, i.e., (3, and 6e, the uncertainties in the 
astrometrically determined values of the angular Einstein ring radius, A(6 *e/6*e,o), will 
depend on the lensing parameters. The values of the lensing parameters ^e and 6e are 
related to the physical parameters of lenses, Doi, v, and M, and thus the uncertainties 
will also depend on the lens parameters. In addition, due to the differences in the lens 
parameters for different types of lensing events, the expected values of A(6 *e/6'e,o) foi’ EMC 
events will be different from those of Galactic bulge events. As a result, the expected 
uncertainties of 6e will have a wide range. 

Boden, Shao, & Van Buren (1998) performed a detailed study of the performance of 
the astrometric observations of gravitational microlensing events in determining various 
observables, including Be- From this analysis, they demonstrated the usefulness of 
astrometric observations of lensing events in determining the lens parameters. In their 
analysis, however, Boden et al. (1998) considered only LMC events, while the majority 
of events are detected toward the Galactic bulge. In addition, their estimate of the 
uncertainties is based on the very limited number of test cases they considered. For 
example, all of the events they considered had a hxed lens location of Doi = 8 kpc, and 
their only values for the Einstein timescale and impact parameter were tE = 0.1 and 0.2 yr 
and f3 = 0.4 and 0.8, respectively. They also did not take into account the restriction the 
source star brightness can impose on astrometric observations. Therefore, one cannot make 
general conclusions about the usefulness of astrometrically monitoring lensing events based 
on their analysis alone. 

In this paper, we investigate the applicability of the astrometric method to the 
determination of the lens parameters for events toward both the LMG and the Galactic 
bulge. For this analysis, we investigate the dependency of A(6 'e/6'e,o) on the lens parameters 
by testing various types of events. In addition, by computing A(6 *e/6'e,o) for events with the 
most probable lensing parameters for a given lens mass under the standard models of the 
Galactic matter density and velocity distributions, we determine the expected distribution 
of the uncertainties in 9e as a function of a lens mass. From this study, we hnd that the 
values of the angular Einstein ring radius are expected to be measured with uncertainties 
A(6 'e/6'e,o) ^ 10% up to the lens mass of M ~ 0.1 Mq for both Galactic disk-bulge and 
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halo-LMC events with a moderate observational strategy. The uncertainties are relatively 
large for Galactic bulge-bulge self-lensing events, A(6*e/6*e,o) ~ 25% for M ~ 0.1 Mq, but 
they can be substantially reduced by adopting more aggressive observational strategies. We 
also hnd that although astrometric observations can be performed for most photometrically 
detected Galactic bulge events, a signihcant fraction (~ 45%) of LMG events cannot be 
astrometrically observed due to the faintness of their source stars. 


2. Dependency of Lens Parameters 


To show the applicability of the astrometric method to the determination of 6e, we 
begin our analysis by investigating how the uncertainties in an astrometrically determined 
6e changes with respect to varying lens parameters. We determine the uncertainties in 6^ 
by conducting model hts to the astrometric shifts of simulated events with various lensing 
parameters. The result of the £t is obtained by computing i.e.. 
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where Wbs is the number of the astrometric observations and (S0co,x, ^^co,y) and {S9c^x, ^^c,y) 
are the centroid shifts of the simulated and model events, respectively. Following the mission 
specihcations of SIM, we adopt the positional accuracy of the astrometric measurements 
to be as 9 ^ = 0.01 mas ([http: / /huey.ipl.nasa.gov/sim|) 


Since our goal is to see the 


variation of A(0 e/^e,o) with respect to individual lensing parameters, the uncertainties 
are determined under fixed but realistic observational conditions. We will discuss the 
expected uncertainties under other observational conditions in § 4. To determine the 
uncertainties, we assume that the astrometric centroid shifts of events are measured with 
a frequency / = 1 day“^ during tmm A Gbs A ^max, where fmin and fmax are the times 
of the hrst and last astrometric measurements. Since the astrometric observations will 
be performed for events that are photometrically detected, we set train = —0.4tE while 
^max = S.OtE- The uncertainties are determined by the 3a level, which is equivalent to 
~ Xmin — 9, where Xmm is best-fitting value. To illustrate the sensitivity 
of this uncertainty on the lens parameters, we have tested a sample event with lensing 
parameters {l3,tE,dE) = (0.5,11.3 days, 0.22 mas). Then the dependency of the uncertainty 
on each lensing parameter is obtained by varying the parameter of interest while holding 
the other lensing parameters constant. 


In Figure 1, we present the uncertainties of the astrometrically determined angular 
Einstein ring radius with respect to various lensing parameters. From the figure, one finds 
that the uncertainties increase significantly with decreasing angular Einstein ring radius 
and decreasing Einstein timescale. On the other hand, the dependency of A(0 e/^e,o) on the 
impact parameter is negligible. 
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3. Dependency on Lens Masses 


In the previous section we showed that the uncertainty A(6*e/6*e,o) depends strongly on 
the size of the angular Einstein ring radius and the duration of the Einstein timescale. Since 
both 9e and ^e are related to the mass of the lens, the strong dependency of A(6*e/6'e,o) 
on 6e and ^e implies that the uncertainties for events caused by different populations of 
lenses will be greatly different. Then the naturally arising questions are: “Is the astrometric 
observation of lensing events a guaranteed method of determining 0 e?” and if not, “To 
what extent can one determine 6e with this method?” To answer these questions, one must 
determine the expected value of A(6'e/0e,o) as a function of lens mass. 

We determine the relation between A(0 e/^e,o) and the lens mass by computing the 
uncertainties for events with lensing parameters which are the most probable for a given 
lens mass. Because the values of 6 *e and ^e also depend on the lens parameters v and Doi, 
there is no one-to-one correspondence between the lens mass and the lensing parameters. 
However, with models of the Galactic matter density and velocity distributions, one can 
statistically determine the most probable values of the lensing parameters corresponding to 
individual lens masses. 


With models for the Galactic matter density and velocity distributions, the distributions 
of tE and 9e for a given lens mass M are obtained by 
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for the angular Einstein ring radius. Here p{Doi) and p{Dos) are the density distributions of 
lens and source stars, 5 is the delta function, Vy and Vz are the components of the transverse 
velocity v, and f{vy,Vz) is their distribution. In the above equations, the factors 7rr| and 
V are included because events with larger cross-sections and higher transverse speeds are 
more likely to occur. For Galactic bulge events (disk-bulge plus bulge-bulge self-lensing 
events), Vy and Vz represent the velocity components which are normal and azimuthal to 
the Galactic plane, respectively. On the other hand, since we adopt a non-rotating isotropic 
velocity model for LMG events (halo-LMG events, see Table 1), these variables represent 
two velocity components of arbitrary direction in the plane normal to the line of sight 
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toward the LMC. We assume that the individual components of the transverse velocities 
have Gaussian distributions of the form 


f{vi) oc exp 



2 < 


i e {y,z}. 


(3.3) 


The values of the mean, Uj, and the dispersion, a^., of the velocity distributions, which are 
listed in Table 1, are adopted from the models of Han & Gould (1995) for the Galactic 
bulge events and from Han & Gould (1996) for the LMG events. For the matter density 
distribntions of the individnal Galactic components, we adopt an axisymmetric model 
(Kent 1992) for the Galactic bulge, a double-exponential disk model (Bahcall 1986) for the 
Galactic disk, and an isothermal sphere model with a core radius (Bahcall, Schmidt, & 
Soneira 1983) for the Galactic halo. The adopted models of the matter density distribntions 
are listed in Table 2. The adopted distances to the Galactic center and to the LMG are 
8.0 kpc and 55 kpc, respectively. 

In Figure 2, we present the expected distribution of /(6*e) and /(^e) for various 
types of events and different values of lens mass. From the hgure, one can see that the 
expected distribntion for a given type of event covers a wide range, considering the variety 
of lens populations and corresponding lens masses. In addition, even with a hxed lens 
mass, the distribntions for different types of events are substantially different from each 
other. Therefore, the analysis of A(6'e/^e,o) based on limited sets of lensing parameters for 
only a single type of lensing events will lead to erroneous conclusions about the general 
applicability of the astrometric method to the determination of 9^. 

Once the most probable values of 9-^ and for a given lens mass are obtained from 
their distributions, the expected valnes of A(6 *e/6'e) are determined as before, by carrying 
out a £t to the astrometric centroid shifts of the event with the most probable lensing 
parameters. Since the dependency of the uncertainties on the impact parameter is not 
important, we assume (3 = 0.5 for all simulated events. The astrometric centroid shifts 
of the events are assumed to be measured under moderate observational conditions with 
a frequency / = 0.5 day“^ dnring —O.d^E < Gbs < 3.0tE and with a positional accnracy 
of = 0.01 mas. Then, the relation between A(6*e/6*e) and lens mass is obtained by 
repeating the nncertainty determinations for different lens masses. 

In Figure 3, we present the resulting relation between the lens mass and the expected 
uncertainty in the astrometrically determined values of 9e- We hnd that the angular Einstein 
ring radii are expected to be measured with uncertainties A(6*e/6*e,o) ^ 10% np to the lens 
mass of M ~ 0.1 Mq for both Galactic disk-bulge and halo-LMG events. For Galactic 
bulge-bulge self-lensing events, the uncertainty is relatively large: A(6*e/6*e,o) ~ 25% for 
M ~ 0.1 Mq. 



4. Improved Astrometric Observational Strategies 


Up to now, we have investigated the uncertainties A(6'e/6*e,o) for a hxed observational 
strategy. In this section, we investigate the dependency of the uncertainties on observational 
strategies to hnd the optimal astrometric observational strategies for better determinations 
of Be- The accuracy of the astrometric determination of 6e can be broadly improved in three 
ways. First of all, the uncertainties will be decreased if the astrometric positional accuracy 
can be improved. However, since the positional accuracy is restricted by the instrumentation 
and not by the observational strategy, we do not consider this scenario further. Secondly, 
the accuracy can be improved by increasing the duration of the astrometric observations. 
Finally, increasing the frequency of astrometric observation will also help to improve the 
accuracy in determining Be- 

In Figure 4, we present the uncertainties as functions of the frequency (upper 
panel) and the duration (lower panel) of the astrometric measurements of 6Bc. The 
uncertainties are obtained for an example event with a set of lensing parameters 
{P,tE,BE) = (0.5,11.3 days, 0.22 mas) under various observational conditions. From the 
hgure, one hnds that the accuracy improves signihcantly with increasing observational 
frequencies. Longer observations of the event also contribute to decrease the uncertainty of 
Be determination. However, observations longer than fmax ~ S^e do not help to a signihcant 
reduction of the uncertainty. 

As we know that measuring SBc with increased frequency will help to improve the 
accuracy in determining Be, we re-determined the relation between A(6*e/6*e,o) and the 
lens mass for a new observational strategy with / = 3 day“^. Since we already tested the 
uncertainties for a long enough observational duration of tmax = 3tE, we only increased 
the observational frequency for this test. In Figure 2, we present the newly determined 
^{Be/Be,o)-M relation (represented by a thick solid line) for the Galactic bulge-bulge 
self-lensing events, which have the largest uncertainties among the various types of events. 
One hnds that the accuracy of the Be determination signihcantly improves - the uncertainty 
decreases by approximately one-half over the entire range of lens masses. 


5. Restriction by Source Star Brightness 

The applicability of the astrometric method in determining lens parameters is 
additionally restricted by source star brightness because the astrometric observations will 
only be possible for events with bright source stars. For both the Galactic bulge and 
LMG helds, the photometric monitoring of source stars is restricted by crowding. Gurrent 
microlensing experiments toward the Galactic bulge reach the photometric detection limit 
when the stellar density of the helds arrives at ~ 10"^ stars deg” (G. Alcock 1997, private 
communication). Based on the model luminosity function of Han (1997), this corresponds 
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to ~ 20.5, which approximately coincides with the astrometric detection limit of 
SIM {y ~ 20, [http:/ /huey.,ipl.nasa.gov/sim|) . Therefore, the astrometric observations of 
microlensing events can be achieved for most of the photometrically detected Galactic 
bulge events. By using the i?-band luminosity function provided by Gould (1998), we 
also determine the fraction of LMG events with source stars bright enough for astrometric 
observation. With the photometric determination limit of G ~ 21 (Alcock et ah 1997b), 
which corresponds to i? ~ 21.1, we hnd that nearly half (~ 45%) of LMG source stars are 
fainter than the astrometric detection limit {R ~ 20.25). This implies that for LMG events 
the major restriction in the general applicability of the astrometric 0 e determination comes 
from the faintness of the source stars and not from the large uncertainties in the determined 
values of 6 *e. 


6. Conclusion 

In this paper, we have considered the effectiveness of astrometric monitoring of 
gravitational microlensing events in determining the lens parameters. The main results 
from this analysis can be summarized as follows: 

1. The uncertainties of the astrometrically determined angular Einstein ring radii are 
strongly dependent on the values of 6 *e and for the events. However, the dependency 
of A(0 e/^e,o) on fhe impact parameter is not important. 

2. The distributions of ^e and ^e for events caused by various populations of lenses 
cover wide ranges. Even with a hxed lens mass, the distributions for different types 
of lensing events are substantially different from each other. Therefore, the analysis 
of A(6 *e/6'e,o) based on limited sets of lensing parameters for only a single type of 
lensing events will lead to erroneous conclusions about the general applicability of the 
astrometric method in determining lens parameters. 

3. Measurements of 59^ with increased frequencies during longer times of observations 
help to improve the accuracy of 0 e determination. However, measurements longer 
than tmax ~ S^E do not contribute to a signihcant reduction in A(6 'e/6*e,o)- 

4. With a moderate observational strategy, the value of 9-^ can be determined with an 
uncertainty A(6 *e/6'e,o) ^ 10% up to a lens mass of M ~ 0.1 Mq for the Galactic 
disk-bulge and halo-LMG events. The uncertainties for the Galactic bulge-bulge 
self-lensing events is relatively large with A(6 'e/6'e,o) ~ 25% for M ~ 0.1 Mq. 
However, the uncertainty can be reduced substantially by adopting a more aggressive 
observational strategy. 

5. Astrometric observations will be possible for most photometrically detected Galactic 
bulge events. However, nearly half of LMG events will not be astrometrically 
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observable due to the faintness of their source stars. Therefore, for LMC events the 
major restriction in the general applicability of the astrometric determination of 9-^ 
comes from the source star brightness and not from the large uncertainties in the 
determined value of 6 *e. 


We would like to thank P. Martini for careful reading the manuscript. We also would 
like to acknowledge the support of the Korea Research Foundation for the 1999 program 
year. 



-11 - 

REFERENCES 

Alard, C., & Guibert, J. 1997, A&A, 326, 1 
Alcock, C., et al. 1997a, ApJ, 479, 119 
Alcock, C., et al. 1997b, ApJ, 486, 697 

Allen, R., Shao, M., & Peterson, D. 1998, Proc. SPIE, 2871, 504 
Ansari, R., et al. 1996, A&A, 314, 94 
Bahcall, J. N. 1986, ARA&A, 24, 577 

Bahcall, J. N., Schmidt, M., & Soneira, R. M. 1983, ApJ, 265, 730 
Boden, A. F., Shao, M., & Van Buren, D. 1998, ApJ, 502, 538 
Gould, A. 1994, ApJ, 421, L71 

Gould, A. 1998, ApJ, submitted (preprint bstro-ph/9807350]) 

Han, G. 1997, ApJ, 484, 555 

Han, G., & Gould, A. 1995, ApJ, 447, 53 

Han, G., & Gould, A. 1996, ApJ, 473, 233 

Hpg, E., Novikov, I. D., & Polnarev, A. G. 1995, A&A, 294, 287 

Jeong, Y., Han, G., & Park, S.-H. 1998, ApJ, submitted (preprint |astro-ph/9807165|) 

Kent, S. M. 1992, ApJ, 387, 181 

Loeb, A., & Sasselov, D. 1995, ApJ, 449, L33 

Maoz, D., & Gould, A. 1994, ApJ, 425, L67 

Nemiroff, R. J., & Wickramasinghe, W. A. D. T. 1994, ApJ, 424, L21 

Paczyhski, B. 1998, ApJ, 404, L23 

Udalski, A., et al. 1997, Acta. Astron., 47, 169 

Walker, M. A. 1995, ApJ, 453, 37 

Witt, H. J., & Mao, S. 1994, ApJ, 430, 505 


This manuscript was prepared with the AAS lATf^X macros v4.0. 







- 12 - 

TABLE 1 

The Transverse Velocity Distribution Models 


event type 

Vy 

(km s ^) 

(km s“^) 

Vz 

(km s~^) 

(km s~^) 

bulge-bulge 

-220(1 -T>) 

[100^(1 + T))2]V^ 

0 

[1002(1 + T))2]V^ 

disk-bulge 

220D 

[302 + 1002T)2]V^ 

0 

[202 + 1002 T ) 2 ] 1/2 

halo-LMG 

0 

250/^2 

0 

250/^2 


NOTE.— The transverse velocity distribution models for different types of 
microlensing events. The distributions are assumed to be Gaussian in form, 
f{vi) = exp —{vi — Viy/2al. , i G {y,z}, and the values of the mean and the 
standard deviation of the distributions are listed. Here D = Doi/Dos- 


TABLE 2 

The Model Matter Density Distributions 

Galactic distribution 

components (Mq pc“^) 

disk p{R, z) = 0.06 exp {— [{R — Rq)/3500 + z/32h\] 

bulge p{s) = 1.04 X 10®(s/0.482)“^-®® (for s < 938 pc) 

p{s) = 3.53if’o(s/667) (for s > 938 pc) 
halo p{r) = 7.9 x 10“^(r^ + RQ)/{rl + r^) 


NOTE.— In the Galactic bulge model, + (2:/0.61)"^, R = {x^ + where 

X and represent the axes directed along the line of sight and toward the Galactic pole. 
The notation Kq represents a modihed Bessel function. The adopted values of the solar 
Galactocentric distance and the core radius of the halo is i?o = 8 kpc and = 2 kpc, 
respectively. 
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Figure 1: The dependency of the uncertainties of the astrometrically determined 0 e on various 
lensing parameters. The astrometric centroid shifts of the event are assumed to be measured 
with a frequency / = 1 day“^ during —0.4tE < tobs < 3.0tE and with a positional accuracy 
C(56»c = 0-01 mas. To investigate the dependencies, we test an example event with a set of lensing 
parameters of (/3, tEj^s) = (0.5,11.3 days, 0.22 mas). The dependency of the uncertainties on 
each lensing parameter is obtained by varying the parameter of interest while holding the other 
lens parameters constant. One finds that the uncertainties A(0 e/0£;,o) increase significantly with 
decreasing angular Einstein ring radius and with decreasing Einstein timescale. On the other hand, 
the dependency of A{9e/9e,o) on the impact parameter is negligible. 
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Figure 2: The distributions of the Einstein ring radii and the Einstein timescales for various types of 
events and different values of lens mass. Erom the figure, one finds that the expected distribution 
for a given type of event covers a wide range, considering the variety of lens populations and 
corresponding lens masses. In addition, even with a fixed lens mass, the distributions for different 
types of events are substantially different from each other. Therefore, the analysis of A(0 r/6*e,p) 
based on limited sets of lensing parameters for only a single population of lensing events will 
lead to erroneous conclusions about the general applicability of the astrometric method to the 
determination of 0 e- 
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Figure 3: The relation between the lens mass and the expected uncertainties in the astrometrically 
determined values of angular Einstein ring radii for various types of events. The uncertainties are 
determined by carrying out fits to the astrometric centroid shifts of the events with lensing 
parameters 6e and which are the most probable for a given lens mass under the models of 
Galactic matter density and velocity distributions. Since the dependency of the uncertainty on 
the impact parameters is not important, we assume f3 = 0.5 for all events. For the moderate 
observational conditions, the astrometric centroid shifts of the events are assumed to be measured 
with / = 0.5 day“^ during —0.4tE < tobs < S.Ots and a positional accuracy of ase^ = 0.01 mas. The 
A(0R/0E,p)-Tf relation for the Galactic bulge events is re-determined under a new observational 
strategy with / = 3 day“^ and it is represented by a thick solid line. 
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Figure 4: The dependency of the astrometrically determined 0 e on the observational frequencies 
(upper panel) and the duration of measurements (lower panel). The uncertainties are determined 
for an example event with a set of leasing parameter of (/3, tEj^s) = (0.5,11.3 days, 0.22 mas) 
under different frequencies and durations of astrometric measurements of the centroid shifts. 
From the figure, one finds that the accuracy of 0 e determination improves significantly with 
increasing observational frequencies. Longer observations of the event also contribute to decrease 
the uncertainty of 9e determination. However, measurements of 56c longer than tmax ~ 3tE do not 
help to a significant reduction of the uncertainty. 





